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ABSTRACT 

Using the Infrared Spectrograph on board the Spitzer Space Telescope, we present low-resolution 
(64 < X/SX < 124), mid- infrared (20-38 /im) spectra of 23 high-redshift ULIRGs detected in 
the Bootes field of the NOAO Deep Wide-Field Survey. All of the sources were selected to have 
1) / 1/ (24/im) > 0.5 mJy; 2) R — [24] > 14 Vega mag; and 3) a prominent rest-frame 1.6 /im stellar 
photospheric feature redshifted into Spitzer 's 3-8 /im IRAC bands. Of these, 20 show emission from 
polycyclic aromatic hydrocarbons (PAHs), usually interpreted as signatures of star formation. The 
PAH features indicate redshifts in the range 1.5 < z < 3.0, with a mean of (z) = 1.96 and a dispersion 
of 0.30. Based on local templates, these sources have extremely large infrared luminosities, compara- 
ble to that of submillimeter galaxies. Our results confirm previous indications that the rest-frame 1.6 
/im stellar bump can be efficiently used to select highly obscured starforming galaxies at z ~ 2, and 
that the fraction of starburst-dominated ULIRGs increases to faint 24 /im flux densities. Using local 
templates, we find that the observed narrow redshift distribution is due to the fact that the 24 /im 
detectability of PAH-rich sources peaks sharply at z = 1.9. We can analogously explain the broader 
redshift distribution of Spitzer-detected AGN-dominatcd ULIRGs based on the shapes of their SEDs. 
Finally, we conclude that z ~ 2 sources with a detectable 1.6 /im stellar opacity feature lack sufficient 
AGN emission to veil the 7.7 /im PAH band. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: formation — galaxies: starburst 
— infrared: galaxies 
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1. INTRODUCTION 

The launc h of the Spitzer Space Telescope 
(| Werner et al.l l2004af ) has allowed the identification 
and study of significant populations of distant, infrared- 
bright galaxies. The most extreme of these are 
exceptionally faint in the optical but readily detected 
in 24 lira surve ys carried out with Spi t zer /M IPS 
(|Rieke et alJl2004| ). For example, iDev et al.l (|2008f l se- 
lected a population of Dust-Obscured Galaxies (DOGs) 
from the Bootes field of the NOAO D eep Wide-Field 
Survey fNDWFS: |ja"nnuzi fc Devlll999l) via the criteria 
R - [24] > 14 Vega mag and f,;(24/im) > 0.3 mJy (see 
also iHouck et al.l l2005t iFiore et all [2008;). Drawing on 
Svitzer /IRS (jHouck et al.ll2004h redshifts (this work and 
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IHouck et all [2005h near-infrared groundb ased spectra 
(iBrand et alj 2*007al). a nd optical spectra (jDesai et al.1 



120081 ). IDev et all (|2008h find that DOGs have a broad 
redshift distribution centered at (z) = 1.99. These 
redshifts imply enormous luminosities, similar to and 
even exceeding those of local ULIRGs and distant 
submillimeter galaxies. While DOGs are rare («2600 
are fo und over the ~9 square degrees of the Bootes 
field), IDev et all (|2008l ) estimate that they contribute 
a quarter of the infrared luminosi ty density at z = 2. 
In addition, their spac e densities (lDev et al.ll2008f ) and 
clustering properties (|Brodwin et alj l2008f ) have led 
to the suggestion that DOGs may be evolutionarily 
related to both coeval SMGs and the most massive local 
galaxies. 

However, the properties of DOGs must be studied more 
thoroughly before any possible relationships to other 
populations can be firmly established. One of the most 
basic unresolved issues is whether the enormous lumi- 
nosities of DOGs are powered predominantly by star 
formation or A GN activity. Th e optical through mid- 
infrared SEDs (|Dev et al.l l2008f) of the DOGs that are 
brightest at 24 /im tend to resemble power-laws, as ex- 
pected for AGN-dominated sources. In contrast, the 
fainter DOGs include a larger fraction of sources featur- 
ing a bump at rest-frame 1.6 /im. This bump is charac- 
teristic of old stellar populations. Its detectability indi- 
cates limited AGN activity because an AGN would result 
in extra flux at rest-frame «2.5 /im, thereby masking the 
bump. 

Follow-up data for the bright (/„(24/im) > 0.75 mJy) 
power-law DOGs support the interpretation that they 
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are AGN-dominatcd. IRS spectroscopy reveals absorbed 
power-laws in the mid-infrared, as exp ected for AGN 
(jHouck et aLll2005l ; IWeedman et al.ll2006D . Near-infrared 
spectroscopy of a small sample of bright pow er-law DOGs 
indic ates that they host powerful AGN (jBrand et alJ 
[20061 ). Similarly, the far-infrared SEDs of bright DOGs 
are similar to the A GN-dominated local ULIRG Mrk 231 
(|Tvler et al.ll2009f) . The morphologies of bright power- 
law DOGs have also been examined dMelbourne et aLl 
120091 12008t iBussmann et~aTl l2009t iDasvra et al.ll2008fl . 
and they tend to be more compact than their fainter 
counterparts, as expected for AGN versus starforming 
regions. 

Because they are relatively more difficult to observe, 
there has been less follow-up of the faint (/i / (24^m) < 
0.75 mjy) bump DOG population. A stacking analy- 
sis shows that their far-infrared properties ar e consistent 
with star-formation, and similar to SMGs (|Pope et al.l 
120081) . However, X -ray stacking analyses have pr ovided 
mixed results (|Fiore et al.ll2008l : iPope et al.ll2008fl. Not 
many IRS spectra of faint bump DOGs exist. IPope et alJ 
(2008) found that 12 out of 70 DOGs in the GOODs 
field have serendipitous IRS spectra. The sources lie in 
the range 0.2 < / l/ (24/xm)/mJy < 1.5. Of these, half 
show PAH features with equivalent widths suggesting 
that they are dominated by star formation in the mid- 
infrared. The PAH-rich sources tend to be the fainter 
ones (/^(24/xm) < 0.7 mjy) and have IRAC SEDs that 
deviate from a power law. In addition, IRS spectra exist 
for infrared bright galaxies with a range of 24 /xm flux 
densities that display the 1.6 /jm stellar bump but do 
not necessarily meet the DOG criterion because they are 
too bright in the optic al (|Huang et al.1 2009: F arrah et al.1 
120081: lYan et al.ll2007h . These spectra indicate that the 
presence of the bump is correlated to the presence of 
PAH features in the IRS spectrum, suggesting that bump 
sources have mid-infrared emission dominated by star 
formation. 

Given the potential importance of the DOG population 
in general, and the fact that the DOG population grows 
with decreasing 24 /im flux density, detailed study of 
the faint bump DOGs is necessary. The first step in 
this process is building a statistically significant sample 
of bump DOGs with redshifts for follow-up study (for 
example, with Herschel). Here we present IRS spectra 
for an additional 23 faint bump DOGs. The 20 of these 
for which we were able to determine redshifts make up 
~ 23% of the 86 sou rces used in the redshift distribution 
of iDev et al.l (l200l . 

This paper is organized as follows. In Sj2j we describe 
our observational data and selection critera. In $3] wc 
present our results, namely the redshifts of our targets, 
the composite IRS spectrum, and their bolometric lumi- 
nosities. We discuss these results in [j4] and summarize 
in EJSJ In the following, we use H = 70 km s" 1 Mpc -1 , 
Ct m = 0.3, and fl A = 0.7. 

2. OBSERVATIONAL DATA AND TARGET SELECTION 

We targeted galaxies for mid-infrared IRS spectroscopy 
based on multiwavelength imaging of the 9.3 deg 2 Bootes 
field of the NDWFS. In the following subsections, we 
describe the survey data, our strategy for selecting high- 
rcdshift starforming ULIRGs, the IRS spectroscopy of 
these ULIRGs, and follow-up observations at 70 and 160 



2.1. Multiwavelength Imaging of the Bootes Field of the 
NO AO Deep Wide-Field Survey 

The 9.3 deg 2 Bootes field of the NDWFS 12 has been 
imaged in the B\y, R, I, and K bands down to 5a 
point-source depths of «27.1, 26.1, 25.4, and 19.0 Vega 
mags, respectively. Additional imaging at the J and 
K s bands was obtained for 4.7 deg 2 of the Boote s 
field through the FLAMEX survey (jElston et al.ll2006f) . 
Approximately 8.5 deg 2 of the Bootes field has also 
been mapp ed (PIP 30) wit h the Infrared Array Cam- 
era fIRAC: lFazio et aLll2004f) on board the Spitzer Space 
Telescope (jWerner et al.l l2004a[ ) . The 5a point-source 
depths of the IRAC Shallow Survey are 6.4, 8.8, 51, 
and 50 iiJy at 3.6, 4 .5, 5.6, and 8 /im, respectively 
(jEisenhardt et al .1120041 ). Approximately 8.74 deg 2 of the 
Bootes field has been imaged with t he Multiband Imag - 
ing Photometer for Spitzer (MIPS: iRTeke et al.l [2004h . 
The la point-source depths of the MIPS survey are 0.051, 
5, and 18 mjy at 24, 70, and 160 /im, respectively. Ap- 
proximately 7 deg 2 have been observed at 1.4 GHz with 
the Westerbork Synthesis Radio Telescope (WRST). The 
data arc characterized by a 13" x 27 " beam and a 
1-dflM s limiting sensitivity of 28 fiJy (|de Vries et al.1 
l2002t ). Although both a reduced mosaic and a catalog 
have been made publically available 13 , we performed our 
own photometry on the reduced mosaic to ensure proper 
dcblending for each IRS target. 

2.2. Selection of high-redshift starforming ULIRGs 

Our goal was to select starforming ULIRGs at z w 2 
for follow-up mid-infrared spectroscopy with the IRS. 
Based on previous mid-infrared, near-infrared, and op- 
tical spectroscopy, we have established that a selection 
criterio n of R — [24] > 1 4 Vega mag results in sources at 
z k. 2 (IDev et al.l 120081: iDesai et all 120081: iBrand et all 
I2007U iHouck et all l2005l )~ Given both the sensitiv- 
ity limits of the IRS and our desire to select high- 
luminosity sources, we only considered objects satisfy- 
ing /j,(24/xm) > 0.5 mjy. Finally, we chose sources that 
display a rest-frame 1.6 /mi stellar photospheric feature 
in their IRAC SEDs. A prominent 1.6 fj,m bump sug- 
gests a limited contribution from AGN- heated hot dust, 
implying that the luminosity of the source is primarily 
generated by star formation. In order to choose the best 
candidates, we eliminated sources that were undetected 
in any of the IRAC channels. 

We adopted Arp 220 as a template for identifying 
sources with a prominent 1.6 fim stellar bump. We var- 
ied only the template normalization and redshift in order 
to fit our Arp 220 template to the four IRAC flux den- 
sities of our sample sources. We also fit a power-law to 
these flux densities. If the best power-law model pro- 
vided a better fit than the best Arp 220 model, we re- 
jected the candidate. If the best-fitting Arp 220 model 
corresponded to a photometric redshift less than 1.5, we 
also rejected the source. 

For the remaining sources, we tested the dependence 
of the best-fit photometric redshift on photometry errors 

12 See |http: / /www. noao.e du/noao / noaodeep / for more informa- 
tion regarding the depth and coverage ot the NDWFS. 
3 http:/ /wwww.astron.nl./wow/testcodc.php?survey=5 
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by performing Monte Carlo simulations. For each of 500 
trials, we randomly perturbed the four observed IRAC 
flux densities. The sizes of the perturbations were based 
on Gaussian distributions centered on the observed flux 
densities and with dispersions equal to the reported la 
photometric errors. In this way, we generated an ap- 
proximate photometric redshift probability distribution 
for each source. 

The full SEDs and photometric redshift probability 
distributions for each source were visually inspected. 
Based on these, a total of 23 sources were selected for 
follow-up IRS spectroscopy. Their positions are listed in 
Table [TJ their photometric properties are summarized in 
Tabled and their SEDs are shown in Figure [TJ The first 
fifteen sources were chosen based upon the tightness of 
their redshift probability distributions (favoring sources 
with the highest signal-to-noise photometry). The last 
eight sources had broad photometric redshift probability 
distributions, but were chosen because they had a signif- 
ic ant probab i lity of lying at z > 2.5. 

iDev et al.l (|2008f) find that the DOG selection picks 
out about 50% of the ULIRGs at z w 2. We now ex- 
plore how our present sample of bump DOGs compares to 
other z ~2 galaxies that have been selected for follow-up 
with the IRS. These include sources det ected in Bootes 
(jHouck et alJl2005l; Wcc dman et al.1120061), the First Look 
burvey f FLS: lYan etTall 120051 iSaiina et all l27)07h and 
SWIRE (|Farrah et all 120081 : iLonsdale et al. [ I2008D . The 
Bootes IRS targets cited above satisfy the same R — [24] 
color cut, but are brighter (/ y (24 pm) > 0.75 mJy) 
than the current sample. Although they were selected 
without regard to IRAC SED shape, the bright 24 /im 
flux density cut resulted in the vast majority having 
power-law SEDs through the IRAC bands. The FLS 
sample consists of bright (/„(24/tm) > 0.9 mJy) 24 /im 
sources with less extreme infrared-to-optical flux density 
ratios than we imposed. Also, the FLS sample satisfied 
a cut of \og w (vS2i/vS%) > 0.5 to try to select for star- 
forming, rather than AGN-like, sources. The SWIRE 
sample was selected to have /„(24/im) > 0.5 mJy and 
strong 1.6 (im stellar bumps peaking in IRAC channel 
2 ("bump-2" sources). Figure [2] shows log-^z/S^/^Ss) 
versus log 10 (^5 , 24/ 'vSr) for these different samples. Our 
sample has similar values of log 10 (1/624 /VS^) compared 
to the FLS and SWIRE samples. This is not surpris- 
ing since all three programs were designed to se lect star- 
forming ULIRGs (see also iBrand et all (l2006h) How - 
ever the brighter Bo o tes sa mples of iHouck et al.1 (|2005l ) 
and IWeedman et al.l (|2006h extend to lower values of 
\og 10 (h , S24 : /vSs). This is consistent with their IRS spec- 
tra being primarily AGN-like. 

IRAC color-colo r diagrams have been used by var- 
ious groups (e.g. lLacv et all 120041 ; IStern et all 120051 ; 
iSajina et al.l l2005l ) to determine the origin of the mid- 
infrared luminosity of galaxies. In particular, these 
groups have defined a wedge-shaped region of this di- 
agram that is populated predominantly by AGN. Fig- 
ure [3] shows one version of the IRAC color-color dia- 
gram and the AGN wedge, along with data from the cur- 
rent sample and for several comparison samples. Most 
of the previously-published Bootes sources have AGN- 
like IRAC colors, consistent with their AGN-like IRS 
spectra. The FLS sample includes some sources in the 
AGN wedge, but has a higher fraction of sources that 



are blue in S%. 0/S4. 5. The former tend to have AGN- 
like IRS spectra, and the latter tend to have PAH-rich 
IRS spectra. Our sample also has very blue S$. 0/S4. 5 
colors. It occupies a similar color-color space as the 
"bump-2" SWIRE sources, which have PAH-rich IRS 
spectra. This plot indicates that the ULIRGs with the 
lowest values of Ss.o/Sk.s are the most likely to have mid- 
infrared luminosities dominated by star formation. The 
sources with low values of S^.o/^.s are those in which 
the 8 fxm flux turns over, indicating the presence of a 
1.6 /Ltm stellar bump and a lack of AGN-heated hot dust. 
Conversely, sources with larger values of Ss.o/S^.s have 
an AGN-heated hot dust component which elevates the 
observed 8 /Ltm flux density. These sources are more 
likely to have mid-infrared luminosities dominated by 
AGN activity. Although we did not explicity include 
a S8.0/S4.5 selection criterion, there are no bump DOGs 
with 7^(24) > 0.5 mJy in our 24 /im catalog that have 
SS.0/S4.5 colors bluer than those in our sample. This 
is likely due to our bump selection criterion, where the 
bump is detected in the 4.5 or 5.8 /im band. 

2.3. Spectroscopic Observations and Data Reduction 

Mid-infrared spectra of the 23 candidates described in 
H2.2I were obtained with the IRS on board Spitzer. The 
z « 2 candidates (the first 15 in Table [1]) were observed 
in program 30334 and the higher redshift candidates (the 
last eight in Table [JJ were observed in program 40441. 
We observed each of our targets with the Long-Low 1 
(LL1) module, which provides wavelength coverage be- 
tween 21 — 34 /.im. At the photometric redshifts of our 
targets, this allows us to detect at least the 7.7 fj,m PAH 
emission complex, which is the strongest of the various 
PAH complexes. 

The IRS has two observing modes: staring and map- 
ping. For faint point sources like ours, mapping mode 
is recommended (see the January 2006 Spitzer memo 
entitled Report on ultradeep IRS spectroscopy of faint 
sources). We therefore used mapping mode to take spec- 
tra of each of our targets in 4 different positions along the 
slit. We first acquired a nearby, bright 2MASS star, and 
then offset to the science target. We chose the longest 
available ramp time (120 seconds) to minimize overhead. 
We used a parallel step size of 33 arcseconds for 4 roughly 
equally-spaced steps along the slit. The number of cycles 
per position ranged from 7 to 18, depending on the 24 
^m flux density and the expected redshift of the source. 
The effective exposure times ranged from 853 to 2194 
seconds. 

The raw data were processed by the S13 pipeline at the 
Spitzer Science Center. The pipeline performed ramp 
fitting, dark sky subtraction, droop correction, linearity 
correction, flat fielding, and wavelength calibration. See 
the IRS Data Handbook 14 for more details on these steps. 
The resulting data products are known as the Basic Cal- 
ibrated Data (BCD). In the following we describe addi- 
tional processing of the two-dimensional spectra (resid- 
ual charge removal, rogue pixel interpolation, and sky 
subtraction) and spectral extraction. 

Although the detector is reset prior to each integra- 
tion, a small fraction (1-2%) of the charge persists be- 
tween frames. If no bright sources are observed, then the 

14 http:/ /ssc. spitzer. caltech.edu/irs/dh/ 
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zodiacal background will be the major contributor to the 
residual charge. Over long AORs, the residual charge can 
build to a significant level, and must be removed. We 
subtracted the background in each BCD by subtracting 
off the median of the counts in pixels 25 through 58 (in- 
clusive) of each row. These pixels were chosen because 
they are unaffected by source flux. 

We interpolated over unstable ("rogue") pixels in the 
background-subtracted frames. A mask of known rogue 
pixels is provided by the SSC, and we identified fur- 
ther bad pixels from the data themselves. We used the 
IRSCLEAN program provided by the SSC to find rogue 
pixels in the 2D data, using the default settings. We fur- 
ther searched for pixels with abnormally high variance 
(> 10a) with time. Pixels were also masked if cither the 
background or the noise surpassed a certain threshhold 
(1000 for both). All masked pixels were interpolated over 
using IRSCLEAN. 

Once the 2D images were cleaned of rogue pixels, we 
created sky images from the data. The background for 
each target position was computed by taking the median 
of the cleaned images for all of the other positions. 

The individual reduced frames were coadded to pro- 
duce final 2D spectra at each map position. One dimen- 
sional spectra were extracted at each map position using 
the SPICE software provided by the SSC. The optimal 
extraction option was used. When this option is set, 
an extraction aperture is not set. Instead, each pixel is 
weighted by its position, based on the spatial profile of a 
bright calibration star. After extraction, we are left with 
4 ID spectra per source. These were averaged together 
to produce the final ID spectrum. 

The signal-to-noisc of the resulting spectra is difficult 
to quantify because of the lack of continuum. We esti- 
mate a noise level of 0.1-0.2 mJy by measuring the dis- 
persion in the residuals of each spectrum after subtract- 
ing off the smoothed composite spectrum (see Section 



2.4. Pointed MIPS 70 and 160 fim observations and 
reduction 

Although the multiwavelength survey of the entire 
Bootes field does include Spitzer MIPS 70 and 160 fim 
observations, these are too shallow to provide interesting 
limits on the far infrared SEDs of the DOGs. The 3-cr 
limits are 15 and 54 mJy at 70 and 160 /xm, respectively. 
We therefore obtained deeper, pointed observations of 
these sources (program 30519, PI Le Floc'h). The data 
were reduced u sing version 3.06 of t he MIPS Data Analy- 
sis Tool fDAD lGordon et al.ll2007T ). None of ou r sources 
were detected at either 70 or 160 /xm. In H3.3I we show 
that the 70 fjxn non-detections are consistent with expec- 
tations based on low-rcdshift templates with IRS spec- 
tra similar to these higher-redshift sources. We used the 
following procedure for determining limits at both wave- 
lengths. 

For the 70 /Ltm data, aperture photometry was per- 
formed at the position of the MIPS 24 fxm source, using 
a 16" radius aperture, an 18" radius inner sky annulus, 
and a 39" outer sky annulus. We applied the aperture 
correction tabulated in the MIPS Data Handbook. We 
also measured the dispersion resulting from performing 
similar measurements over empty areas of the 70 micron 
images. The adopted 3er limit (5 mJy) is taken as three 



times this dispersion. 

For the 160 /xm data, aperture photometry was per- 
formed at the position of the MIPS 24 /xm source, using 
a 32" radius aperture, a 64" inner sky annulus, and a 
128" outer sky annulus. We applied the aperture correc- 
tion tabulated in the MIPS Data Handbook. The limit 
was calculated in the same way as for the 70 ^m data. In 
this case, the limit turned out to be 30 mJy, but was very 
uncertain due to the small number of pixels in the 160 
^m images. To be conservative, we adopted the higher 
limit of 40.8 mJy predicted by the Spitzer SENS-PET. 

3. RESULTS 

The IRS spectra of our 23 targets are shown in Figure 
2J Strong PAH features are common in this sample, espe- 
cially when compare d to the brighter ULIRGs observed 
by th e IRS in Bootes (jHouck et al.ll2005l ; IWeedman et all 
12001 . 

In this section we first describe how we determined 
the spectroscopic redshifts of these sources. We then 
use these spectroscopic redshifts to present the composite 
spectrum of the sources with measured redshifts. From 
the composite spectrum, wc measure the characteristic 
7.7 /jm PAH equivalent width. Finally, we estimate the 
bolometric and rest-frame 24 /xm luminosities of these 
sources, so that they may be easily compared with other 
samples, both coeval and at low rcdshift. 

3.1. Redshifts 

Many of the spectra displayed in Figure [4] contain at 
least one PAH emission feature (centered at rest-frame 
6.2, 7.7, or 11.3 /im) from which we can determine a 
rcdshift. Because these PAH features are fairly broad 
and the spectra are of only moderate signal-to-noise, we 
choose to measure the redshift using templates rather 
than fitting the lines individually. We use four templates 
that contain strong PAH features: the unobscured Galac- 
tic reflection nebula NGC 7023, which represen ts a pure 



photo-dissociation region (I Werner et alj 2004b[) : the av- 
erage starburst spectrum of lBrandl et all (120061); a nd the 
prototypical starbursts NGC 7714 (jBrandl et all l2004h 
and M82 (jSturm et al.ll2000f) . Figure shows an exam- 
ple of our redshift determination. After computing the 
best-fit redshift resulting from each template, we adopt 
the redshift associated with the template that provides 
the best formal fit. In 12 cases, this was NGC 7023; 
in four cases this was the average starburst spectrum; 
in three cases this was NGC 7714; and in one case this 
was the M82 template. We were unable to determine 
a spectroscopic redshift for three sources because they 
had no strong features. The adopted redshift for each 
source is listed in Table Q] For a given source, all of the 
templates provided very similar redshifts. The rcdshift 
errors presented in Table [T] are based on the range of red- 
shifts computed for a given source. The average redshift 
of sources in our sample is (z) = 1.96, and the dispersion 
is 0.30. 

Figure [6] presents the rest-frame SEDs, scaled to 
unity at 1.6 /im, of all 20 sources with spectroscopic 
redsh ifts. For compari son, we ov erplot the SEDs o f 
M82 (jSturm et al] 12000]), Arp 220 (jArmus et al.ll2007h . 
Mrk 231 (|Armus et al.ll2007h . The SEDs of the z « 
2 starforming ULIRGs most closely resemble that of 
Arp 220, in terms of their R - [24] colors. The SEDs 
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of our sample appear to be much more strongly red- 
dened than that of M82. Again, this is unsurprising 
given the selection criterion R — [24] > 14 Vega mag. 
Also as expected, the SEDs of the high-redshift ULIRGs 
are significantly colder than that of Mrk 231. The excep- 
tion is Source 21 at z = 3.04. Its observed-frame 24 /tm 
(rest-frame 5.9) flux density is more consistent with the 
Mrk 231 template than that of Arp 220. 

3.2. Composite spectrum 

Because the signal-to-noisc ratio of any individual 
spectrum in Figure [4] is fairly low, we provide three com- 
posite spectra in Figure [7l the median, the mean, and 
the variance-weighted mean. To compute the composite 
spectra, we shifted the wavelengths of each spectrum to 
the rest-frame, based on the redshifts calculated in §3.1| 
sampled the rest-frame wavelengths to a common grid; 
and normalized each spectrum to unity at 7.7 /im. 

The median spectrum shows strong PAH features at 
7.7, 8.6, and 11.3 /im. In the bottom panel of Figure 
El we compare this median sp ectrum to an averag e star- 
burst spectrum compiled by iBrandl et al.l (|2006f ). to a 
me dian cold (S25/S60 < 0.2) ULIRG spectrum compiled 
bv iDesai et a .1 (l2007al) , and to the spectrum of Arp 220 
(|Armus et al. [20071 ). The wavelength range between 9 
and 11 /im can be strongly affected by a broad silicate ab- 
sorption feature centered at 9.7 /im. Unfortunately, the 
limited wavelength range of our spectra prevent an accu- 
rate determination of the unabsorbed continuum level, 
and we are therefore unable to precisely quantify the 
depth of the silicate absorption. However, Figure [7] sug- 
gests that this absorption is comparable to that of local 
ULIRGs with cold SEDs. 

In Figure [8] we compare the median mid- infrared spec- 
trum of our starforming DOGs to those of other high- 
rcdshift galaxies observed with the IRS (obtained via 
private communication). The first comparison sample is 
a subset of the Spit zer-selected ULIRGs detected in th e 
First-Look Survey (|Yan et al.|[200l ISaiina et al.ll2007l ). 
Of 48 targets spanning the redshift range 1 < z < 3, 
25% show strong PAH features. The average strong-PAH 
spectrum is shown in the second panel of Figure [SJ 

The second comparison sample consists of 24 SMGs (23 
detected), with 19 displayin g PAH features, in the red- 
shift range 0.65 < z < 3.2 ( Menendez-Delmestre et all 
|2009|) . Of these, 19 display prominent PAH emission. 
The composite spectrum from this sample is shown in 
the third panel of Figure [8] 

The third comparison sample consists of 13 GOODS-N 
submillimeter gala xies in the redshift range 0.9 < z < 2.6 
(jPope et alJ I2008T ). Of these, 11 have mid-infrared lu- 
minosities dominated by PAH emission. The compos- 
ite PAH-rich spectrum from this sample is shown in the 
fourth panel of Figure [3 

Although Figure [7] provides the visual impression that 
the z « 2 bump DOGs have mid-infrared spectral prop- 
erties intermediate between those of local cold ULIRGs 
and local starburst galaxies, we wish to measure the 
equivalent widths of the PAH emission to facilitate a 
quantitative comparison to other samples. Three PAH 
emission features are clearly seen in Figure including 
those centered at rest-frame 7.7, 8.6, and 11.3 /mi. The 
11.3 /im feature is strongly affected by the rest-frame 
9.7 /<m silicate feature. The other two features are close 



enough in wavelength that they arc often measured to- 
gether. We adopt this procedure and call the result the 
7.7 /tm PAH equivalent width. Because our spectrum 
has limited wavelength coverage, the continuum is quite 
difficult to determine. We the refore adopt an em pirical 
method similar to that used bv lSpoon et a l. (2007). The 
technique is illustrated in Figure [5] The resulting rest- 
frame equivalent width for the 7.7 /tm PAH feature is 
0.5±0.05 /tm for our sample. Given the difficulties of this 
measurement, this is similar to the values of 0.4±0.05 /tm 
that we find for the comparison samples. I t is also simi- 
lar to the values found for local starbursts (jBrandl et alJ 
12001 . 

In addition to constructing a single composite spec- 
trum, we also tried binning the spectra by observed prop- 
erties to explore the variation within our sample. For 
example, a trend with spectral type and luminosity is 
observed in the local univers e (see ^4.2[) . It has a lso been 
seen in the z as 2 sample of ISaiina et al.l (|2007l ). How- 
ever, we find no convincing correlation between spec- 
tral shape and iTj„(24/im) in our sample (see M3.3I for 
derivation of luminosity). In both the low and high 
redshift samples where a correlation is seen, the scat- 
ter is very large, and can only be seen when compar- 
ing the full range of spectral types (AGN-dominated to 
starburst-dominatcd) . Therefore, the lack of a correla- 
tion within our modest sample of primarily starburst- 
dominated sources is consistent with previous findings. 
We also explored spectral variations with the colors plot- 
ted in Figures [2] and [3l and found no convincing trends. 

3.3. Luminosities 



Given the redshifts calculated in £13.11 ((z) = 1.96), we 
can estimate the infrared luminosities of our sources. Our 
estimates are necessarily crude because we currently have 
no direct flux measurements at wavelengths that probe 
the expected peak (« rest-frame 70 /tm) of the spectral 
energy distributions of our sources. Our longest wave- 
length photometric detection comes from Spitzer MIPS 
observations in the 24 /tm band, which corresponds to 
a rest-frame of «8.4 /tm for these sources. Although we 
have 70 and 160 /tm observations (corresponding to rest- 
frames 24 and 55 /tm, respectively), these provide only 
limits, not detections (see £ 12. 4j) . 

Given our lack of direct constraints on the far-infrared 
SEDs of our sources, we extrapolate from the rest-frame 



mid-infrared to determine La 



-iooo/im- 



There arc a va- 



riety of bolomet ric correction factor s calibrated on lo- 
cal sources (e.g. ICaputi et al.l I2007D and several fami- 
lies of templates (e.g. Chary fe Elbazl |2001|; iDale et al 
120051 : ISiebenmorgen fe Krugell l2007t iRieke et all l200" ^ 
that can be used to guide this extrapolation. Unfortu- 
nately, the resulting value of Lg_iooo/tm can vary by a 
factor of 5 depending on which local relation or template 
is use d (pale et al.ll2005l : ICaputi et afll2007l : iDev et all 
l2008f h Thus, it is important to choose local templates 
that are as similar as possible to the sample with un- 
known far-infrared SEDs. 

A first pas s wit h the starburst templates of 
I Chary fc Elbal (|2001[ ) indicates that the bump DOGs 
have ULIRG luminosities. To refine our estimate of 
the bolometric luminosity, we therefore wish to use 
a local ULIRG sample for our extrapolation. The 
Spitzer IRS GTO Team obtained IRS spectra for 
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and luminosities 


Armus et al.1 


2007:lDesai et al.ll2007al: 


iFarrah et alj|2007 


; ISpoon et al. 


I2007D. We use this sam- 



pie to estimate Lg_iooo/im using the following procedure. 
First, we fit each of the 107 GTO ULIRGs to the me- 
dian spectrum shown in Figure We then ranked the 
GTO ULIRGs according to how well they fit the me- 
dian bump spectrum, as measured by the % 2 - We fit 
the first-ranked GTO ULIRG to each of the individual 
bump spectra, yielding infrared luminosities in the range 
2.5 < L 8 _iooo J tim/10 12 L Q < 25. This variation is due 
both to the fact that the bump sources have a variety 
of 24 /im flux densities (from 0.5 to 0.75 mJy) as well 
as a range of rcdshifts (from 1.5 to 3.0). To examine 
the possible error in the value of Ls-ioooum for a given 
bump source, we calculated the range in Lg-iooo/im that 
results from the most highly ranked 10, 25, and 50 lo- 
cal GTO ULIRGs. These result in Ls-iooo/iin ranges 
that are within a factor of 2, 3, and 3 of the best-fit 
Ls-iooo/im- Wc estimate that the individual values of 
Ls-iooo/im are good to a factor of «3. The results are 
listed in Table H 

The radio fluxes and limits presented in Table [2] (see 
also t|2.ip can also be used to constrain L 8 _iooo/im- At 
low redshift, a far- infrared-radio correlation has been ob- 
served to hold for a wide variety of galaxy types over four 
orders of magnitude in luminosity (e . g. Ide Jong et al.1 



s ot magnitude m luminosity (e . g. Ide Jong et al.l 
iHelou et al.l fl985t lCondonlfT99a lYun et al.ll200lt 



1985^ 

Bell 120031 ). Such a correlation has also been observed 



at higher redsh i fts (0 . 6 < z < 3; iGruppioni et al.l 12003b 
Appjeton et al.l 12004 iFraver et all 120061: IKovacs et all 



20061) . To illustrate the systematic uncertainties in using 
the far-infrared-radio correlation to estimate Lg_iooo/im 
for our bump sources, we use two different calibra tions 
of the correlation. The first calibration is that of iBelll 
(|2003f ). and is based on a sample of 164 local galax- 
ies without s i gns o f AGN activity. The second is by 
IKovacs et al.l (|2006l ). and is determined from 15 high- 
redshift (1 < z < 3) SMGs. The high-redshift calibration 
yields smaller values (by a factor of «4) of L 8 _i oo/im for 
a measured radio flux. It is consistent with the findings 
of iMurphv etal] (|200l for 11 SMGs out to z = 2.5. The 
infrared luminosities resulting from each calibration are 
presented in Table [3] Despite the large uncertainties, 
the radio-derived luminosities are consistent with the in- 
terpretation that these sources are ULIRGs. Because a 
significant fraction of our sources have only radio lim- 
its, in the following we adopt the values of L 8 _iooo/im 
derived by extrapolating the IRS spectra of our sources 
using local ULIRG templates (column 2 of Table [3]). 

If all of this luminosity is produced by star formation, 
this would correspon d to star formatio n rates in the range 
400-4000 Mq yr" 1 (jKennicuttl [l998h . In reality, some 
fraction of this infrared luminosity is likely due to AGN 
activity. Unfortunately, given the limited wavelength 
range of our spectrum, we cannot disentangle the frac- 
tion of the mid-infrared luminosity that is contributed 
by star formation versus AGN activity. However, the 
AGN contribution is likely not dominant, since a strong 
power- law component would overwhelm both the 1.6 /im 
stellar bump and the PAH features. Although these es- 
timates are uncertain, they indicate that the sources we 
have identified arc extremely luminous, comparable to 
the coeval submillimeter galaxies. Another possibility is 



that the local templates differ significantly from our high- 
redshift sources in the far-infrared. Direct measurements 
at long wavelengths are needed to explore this possibility. 

In ij2.41 wc described follow-up 70 and 160 /im follow- 
up observations of a subset of the sources presented in 
this paper. All were non-detections, with 3-<r limits of 
5 and 40.8 mJy at 70 and 160 /im, respectively. If our 
sources are similar in luminosity to SMGs, then these 
limits should be consistent with the far-infrared SEDs 
measured fo r SMGs. Using th e SMG SEDs presented in 
Figure 5 of iPope et all (|2008l) . we estimate that SMGs 
with 24 /im flux densities similar to our bump sources 
should have observed 70 /im flux densities of ~2 mJy and 
160 /im flux densities of ~20 mJy. These estimates are 
below our detection limits, and thus our far- infrared non- 
detections are still consistent with SMG-like luminosities. 

Similarly, the wavelength coverage of the available IRS 
spectra allows us to predict the 70 /im flux density that 
would be observed if the local GTO ULIRGs were shifted 
out to 2 w 2. Using the 10 GTO ULIRGs that provide 
the best fit to the median spectrum shown in Figurc[7J wc 
estimate /„(70/mi) =0.5 — 1 mJy. Thus, the observed- 
frame far-infrared SEDs of our bump sources are consis- 
tent with low-redshift ULIRGs. 

4. DISCUSSION 

4.1. Why does our selection of starforming ULIRGs lie 
within such a narrow redshift range? 

Table [TJ shows that our sample of starforming ULIRGs 
has a very narrow range of redshifts: (z) = 1.96±0.3. We 
imposed two selection criteria which impact this distribu- 
tion: (1) our requirement that we see the 1.6 /im stellar 
bump in the IRAC bands; and (2) our requirement that 
/„(24/im) > 0.5 mJy. In this section, we explore how 
these criteria could have resulted in the observed narrow 
redshift distribution. 

Based on local templates, is it plausible that require- 
ment (1) results in the observed narrow redshift range? 
In Figure [31 we showed that requirement (1) translates 
roughly into selecting targets from a region of IRAC 
color-color space defined by log 10 (S' 8 .o//S'4.5) < 0.2 and 
0.0 < log 10 (S , 5.8/S , 3.6) < 0.3. In the top two panels of 
Figure [9l we plot the expected IRAC colors of three lo- 
cal templates as a function of redshift. The three tem- 
plates span a range of properties: Arp 220 is PAH- 
rich ULIRG, 08572+3915 is a deeply obscured AGN- 
dominated ULIRG, and the starburst tem plate is an av- 
erage of lower-luminosity starbursts from iBrandl et al.1 
(|2006l ). Figure [HI shows that Arp 220 meets requirement 
(1) over the redshift range z > 1.75, while 08572+3915 
and local starbursts would never be identified as bump 
sources over the redshift range plotted. We conclude that 
our requirement of seeing the 1.6 /im bump in the IRAC 
bands selects for high redshift sources (z > 1.75), but 
can not by itself explain the narrowness of the resulting 
redshift distribution. 

To investigate how requirement (2) affects the redshift 
range of our sources, we also show the expected 24 /im 
flux density of our three templates as a function of red- 
shift in the bottom panel of Figure [5] All templates 
show a dip at z ~ 1.5. One cause for this dip is that the 
24 /im bandpass is sampling the 9.7 /im silicate absorp- 
tion feature at z w 1.5. In Arp 220 and in the starburst 
template, an additional reason for this dip is that the 24 
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/jm bandpass samples the 8.6 and 11.3 /im PAH emis- 
sion features at both higher and lower redshifts, respec- 
tively. In these two sources, the 7.7 and 8.6 /im PAH 
emission features passing through the 24 /im passband 
cause a peak at z = 1.8 — 1.9. The 24 /mi flux density of 
the AGN-dominated ULIRG 08572+3915 does not fall as 
sharply at z > 1.9 because of strong continuum emission 
longward of rest-frame 8.6 /urn. Also shown in Figure [9] 
is the rcdshift distribution of the bump DOGs with IRS 
redshifts. The redshifts of the bump sources coincide 
with the peak seen in the Arp 220 template. This coin- 
cidence suggests that the narrow redshift distribution of 
the bump sources is due to the fact that they are prefer- 
entially detected over the redshift range where their PAH 
emission features are sampled in the 24 /im bandpass. 

Figure [9] formally allows for the possibility of a sig- 
nificant population of bump DOGs at z > 2.2 that are 
missed by our selection technique because there is no 
PAH feature in the 24 /jm bandpass to boost the 24 /im 
flux density enough to meet our flux cut. However, there 
is a second reason why we do not see bump DOGs at 
z > 2.2. On average, bump sources that satisfy our 24 
(im detection limit at z = 2.2 (z = 2.5) would have to be 
more luminous by a factor of >1.25 (>2) compared to 
those at z = 2. Since the top end of the galaxy luminos- 
ity function typically declines exponentially, this could 
be the reason that we do not see exceptionally luminous 
DOGs powered by star formation at z > 2.2. Thus, both 
the mid-infrared spectral shape of starforming ULIRGs 
and the declining luminosity function at high luminosi- 
ties are responsible for the decline in the redshift dis- 
tribution at z > 2.2 (and the resulting tightness of the 
rcdshift distribution). Which mechanism dominates is 
unclear from the data in hand. 

In the bottom panel of Figure [5] we also show the 
rcdshift distribution of the bright (/ t/ (24 / mi) > 0.75 
mJy) power-law DOGs from iHouck etafl (|2005h . For 
1.3 < z < 2.9, this distribution tracks the profile of 
08572+3915. At lower redshifts, sources become brighter 
in the i?-band, and fail to make our R — [24] DOG cut. 
This suggests that the redshift distribution of power-law 
DOGs is determined in large part by the SED shape 
of AGN-dominated ULIRGs. The fact that power-law 
DOGs are found at greater redshifts than bump DOGs 
suggests that the luminosity function of power-law DOGs 
may be shifted to higher luminosities compared to bump 
DOGs. This is consistent with the observation that the 
most luminous ULIRGs in the local universe tend to 
be AGN-dominated (e.g . IVeilleux et all 19991 : iTran et all 
l200lt iDesai eT"aTll2007al) . 

The presence of a bump in the IRAC SEDs allowed 
the selection of high-rcdshift (z > 1) sources, but the 
additional requirement of a 24 /im detection tightly con- 
strained the redshift distributi on to z ~ 2, for the rea sons 
described above. For example. iBrodwin et all (|2006f l use 
the optical, near-infrared, and IRAC data available in the 
Bootes field to compute photometric redshifts for a sam- 
ple selected at 4.5 /im. They do not use the 24 /im flux 
densities. Using a hybrid method involving both neural 
net and template-fitting techniques (but no ULIRG tem- 
plates), calibrated using over 15,000 spectra, they achieve 
an accuracy of a = 0.06(1 + z) for 95% of galaxies at 
< z < 1.5 and a = 0.12(1 + z) for 95% of AGN at 
< z < 3. However, the photometric redshifts computed 



for the sample presented here range from 1 < z p hot < 3.5, 
with an accuracy of a = 0.8. This shows that photomet- 
ric redshifts that achieve high accuracy for the majority 
of objects do not perform nearly as well on these rare, 
extreme DOGs. 

4.2. Comparison to Local ULIRGs 

We intitially used Arp 220 as a template for selecting 
bump DOGs. Figurc[6]illustrates that the SEDs of bump 
DOGs arc consistent with that of Arp 220. However, 
Figure [9] shows that neither Arp 220 nor 08572+3915, if 
placed at z = 2, would be selected in this study. The 
reason for this is that they would be too faint at 24 //m, 
by factors of 25 and 2, respectively. Even if we scaled 
Arp 220 to match the observed 4.5 /im flux densities of 
our bump DOGs at z = 2, it would not make our 24 /im 
cut. In contrast, a scaled Mrk 231 at z = 2 would. This 
implies that PAH features are not necessary to boost the 
observed-frame 24 /jm flux density of a given source into 
our sample; continuum is sufficient. The fact that we 
see very few continuum-dominated sources in our bump 
sample implies that z w 2 galaxies with only a small 
AGN contribution at observed 8 /im also have a limited 
AGN contribution at observed-frame 24 /im. 

Figure [10] shows the rest-frame equivalent widths of 
the 7.7 and 11.3 /im PAH features as a function of rest- 
frame 24 /im luminosity for the bump sam ple, the z 2 
PAH-r ich Spitzer -selected FLS sample from lSaiina et al" 
J2007D, and for the wlOO local ULIRGs from IDesai et al. 
( 2007al) . This plot illustrates that the population of 
high-redshift starbursts that Spitzer is detecting includes 
galaxies that are uncommonly luminous at rest-frame 24 
/im compared to ULIRGs in the local universe. 

4.3. Why do bump sources have strong PAHs? 

The existing IRS spectra of ULIRGs at z w 2 indi- 
cate that prominent rest frame 1.6 /im stellar photo- 
spheric features are accompanied by large PAH equiv- 
alent widths. The corollary also holds: 24 /im sources 
with power-law rest frame near-infrared SEDs have small 
PAH equivalent widths. 

This trend can be understood within a commonly ac- 
cepted framework, wherein the rest-frame near-infrared 
and mid-infrared spectrum of a galaxy is a combination 
of th ree basic componen ts: stars, the ISM, and an AGN 
(e.g. iHuang et alJ l2007). In galaxies selected by their 
R — [24] colors, the SEDs of these components are also 
strongly affected by dust attenuation. Old stars have an 
SED resembling a blackbody that peaks at 1.6 /an in 
the rest frame and falls off steeply towards longer wave- 
lengths. In contrast, dust-obscured young stars have an 
infrared SED which is characterized by mid-infrared PAH 
features with high equivalent widths and a broadband 
SED which rises steeply beyond 6 /im. Finally, the AGN 
contributes a power-law SED through the near- and mid- 
infared. 

Within this framework, galaxies with an observable 1.6 
/zm stellar bump have a limited AGN contribution. In 
contrast, galaxies without an observable 1.6 fjxa. bump 
have an AGN contribution which is large enough to dis- 
guise the underlying stellar emission by providing addi- 
tional flux on its long- wavelength side (rest frame «2.5 
/tm). If an AGN is emitting significantly at such short 
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wavelengths, and if it has a power-law SED, then the ex- 
pected flux at rest frame wavelengths of «8 /im would be 
even more significant. A strong AGN-heated dust com- 
ponent at 8 /im would increase the continuum under the 
PAH features. The result is a low equivalent width for 
PAH features in sources with weak 1.6 /jm stellar bumps. 
Similarly, a source which does not have enough AGN- 
heated dust to swamp out the 1.6 /im stellar bump will 
not have enough AGN-heated dust to depress the equiv- 
alent width of the mid-infrared PAH features. In this 
scenario, the PAH features are diluted by AGN-heated 
smooth dust continuum. The degree of the dilution de- 
pends on the luminosity and extinction of the AGN rel- 
ative to that of the star burst component. iTeplitz et al.l 
((2007t) mix the SEDs of a local starburst (NGC 7714) 
and a pure QSO (PG 0804+761) to quantitatively ex- 
plore how an AGN contribution affects both the 1.6 /im 
stellar bump and the PAH features. They find that a 
weak AGN contributing 10% of the total 1-1000 /im lu- 
minosity of a composite source can hide the 1.6 /jm bump 
and decrease the PAH EQW by over 40%. 

While it is generally agreed upon that weak PAH fea- 
tures indicate a strong AGN component, PAH dilution 
is not the only scenario that has been invoked as an ex- 
planation. It is also possible that PAH carriers are de- 
stroyed in the harsh radiation field s produced by AGN 
(|Aitken fc Roche) [l98l IVoitl IT992T) . With the data in 
hand, we cannot distinguish between the dilution and 
destruction scenarios within the z m 2 bump sources. 

There is some evidence that PAH dilution occurs in lo- 
cal ULIRGs. Among local cold ULIRGs, the median 11.3 
/urn PAH equivalent width (EW) is 80% of that found in 
starburst ga laxies, while the m edian 6.2 /im PAH EW 
is only 50% (jDesai et al.ll2007bl ). In the destruction sce- 
nario, the EWs of both PAH features should be reduced. 
However, the observation that the 11.3 /im PAH EW is 
more similar to the starburst value finds a natural ex- 
planation in the dilution scenario. Namely, the emission 
from the hot dust responsible for the dilution at 6.2 /mi 
is highly cxtincted at 11.3 /im, due to the strong, broad 
absorption feature centered at 9.7 /im. Whether or not 
these results based on lower redshift ULIRGs apply to 
higher redshift samples remains an open question. 

5. CONCLUSIONS 

We have tested a technique for using mid-infrared pho- 
tometry to select strongly starforming and highly ob- 
scured ULIRGs at z w 2. In a previous work, we demon- 
strated that the simple criteria R — [24] > 14 Vega mag 
and / 1/ (24/im) > 0. 3 mJy select for about half of the 
ULIRGs at z = 2 (jDev et alj [20081 ). We refer to the 
resulting sample as Dust-Obscured Galaxies, or DOGs. 
Initial studies of the DOGs with the brightest 24 /tm 
flux densities (/ 1/ (24/im) > 0.75 mJy) revealed AGN- 
dominated sources. In this paper, we attempted to select 
the most vigorously star-forming DOGs. Our technique 
was to search for DOGs whose IRAC (3-8 /im) SEDs 
show evidence for the redshifted 1.6 /im stellar photo- 
spheric bump characteristic of old stars. Sources with 
strong AGN components should have enough hot dust 
to swamp this feature. Thus, sources with clear 1.6 /tm 
bumps should have a limited AGN contribution to their 
mid- infrared luminosities. These sources are likely to be 
dominated by star formation. In addition to selecting for 



starforming galaxies, this technique also selects for galax- 
ies at z w 2, since the 1.6 /im bump must be shifted into 
IRAC channels 2 or 3 (4.5 or 5.8 /im, respectively) in 
order to be easily detected. To test this selection, we 
obtained mid-infrared IRS spectroscopy of 23 candidate 
star-forming galaxies. Based on these data, we come to 
the following conclusions: 

1. Of 23 targets, 20 have mid- infrared spectra display- 
ing PAH emission lines indicating that 1) these 
sources lie at (z) = 1.96 with a dispersion of 
0.3; and 2) these sources are strongly starform- 
ing. These results indicate that even for the most 
obscured galaxies, the rest-frame 1.6 /im stellar 
photospheric bump can be used to efficiently se- 
lect starforming galaxies at z « 2. This red- 
shift distribution implies that the space density 
of bump DOGs (R - [24] > 14 Vega mag, 0.0 < 

log 10 (5 5 .8/S3. 6 ) < 0.3, -0.3 < log 10 (S 8 . /S 4 .5) < 
0.2) is «7.5 x 10 5 Mpc 3 , consistent with the 
space densit y of SMGs in GO ODS-N (private com- 
munication; (Wall et al.ll2008l ). 

2. We have confirmed that DOGs with faint 24 /tm 
flux densities are more likely to b e star forming than 
brighter DOGs. iBrand et all (|2006f) used X-ray 
data to show that this is true for the gener al popu- 
lation of 24 /tm sources. iDev et all (j2008| ) showed 
that the fraction of DOGs with AGN-like (power- 
law) IRAC SEDs smoothly decreases with decreas- 
ing 24 /tm flux density. In this paper, we verify that 
DOGs which show evidence of the 1.6 /tm stellar 
opacity feature in their mid-infrared photometry 
also show strong PAH emission, and appear to be 
powered principally by star formation. 

3. Local templates predict that our sample of star- 
forming DOGs have infrared luminosities in the 
range L 8 _i 00 o/mi = (6-50) xlO 12 L . If all of 
this luminosity is due to star formation, the star 
formation rates among our sample would be 1000- 
8500 M Q yr" 1 . AGN may be contributing to the 
total infrared luminosity, but the strong PAHs sug- 
gest that star formation dominates at least the mid- 
infrared luminosity. The unphysically large star 
formation rates implied may point to an evolution 
in the template from high to low redshift. This 
is also implied by the fact that starburst galaxies 
as luminous as those recently detectable by Spitzer 
are uncommon in the local universe (|10[). Other 
auth ors have found hints that this may be the case 
fe.g. lDesai et~aDl2007bt [Rigbv et al.ll200l . While 
we have limits on the rest-frame 24 and 55 /tm flux 
densities of these DOGs from MIPS 70 and 160 
/tm data, direct measurements of the far-infrared 
SEDs of high-redshift dusty galaxies are critical to 
determining their contribution to the global star 
formation rate at that epoch. 

4. The redshift distribution of bump DOGs is remark- 
ably narrow because at z = 1.9, the strong 7.7 
/im PAH feature boosts the 24 /tm flux, pushing 
sources with insufficient continuum into our flux- 
limited sample. Analogously, the redshift distri- 
bution of power-law DOGs is dictated by the fact 
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that sources with sufficient continuum to meet the 
24 /xm flux density cut fall out of the sample at 
z w 1.5, when the 9.7 /xm absorption feature is in 
the 24 /xm bandpass. 

5. The dearth of objects with the 1.6 /xm stellar fea- 
ture with AGN-like mid-infrared spectra indicates 
that objects lacking enough AGN emission at 2 /xm 
to hide the stellar bump also lack an AGN powerful 
enough to dilute and/or destroy the PAH features. 
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Fig. 1. — The individual SEDs of the bump DOGs a re sh own as filled points. The red line represents the mid-infrared spectrum of Arp 
220, shifted to the spectroscopic redshift if known (see £|3.1|l and normalized to best fit the IRAC data (red points). The rest-frame 1.6 /^m 
stellar photospheric bump is clearly seen in the Arp 220 spectrum, and is detected in objects at 2 fa 2 as a turnover in the long-wavelength 
IRAC channels. All sources in this paper were selected to exhibit this turnover. No color correction has been applied to the Spitzer 
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spectra: the average starburst template from lBrandl et al.l (1 20061 ) (red); M82 (dark blue); NGC 7714 (green); and NGC 7023 (light blue). 
The redshift associated with the best-fitting template is adopted and the error is taken to be the range of redshifts found for all four 
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0.1 1.0 10.0 

Rest Wavelength (|Lim) 

Fig. 6. — Optical to mid-infrared SEDs of z w 2 bump DOGs. Each SED is arbitrarily scaled to have a 1.6 fim flux density of unity. The 
SEDs of bump DOGs strongly resemble that of Arp 220. As in Figure [T] no color correction has been applied to the Spitzer photometry. 
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Fig. 7. — Top: Median spectrum of 20 z fss 2 bump DOGs. The spectrum is normalized to unity at 7.7 (im. Bottom: Same as top panel, 
but with comparison spectra (also normalized at 7.7 fim) overplotted. The mid-infrared spectrum of bump DOGs appears intermediate 
between those of local starbursts and local cold ULIRGs. 
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Fig. 8. — Illustration of PAH equivalent width measurements and comparison to other high-redshift samples. Panel (a) shows the median 
spectrum of the bump DOGs (black line). To determine the equivalent width of the PAH features in this spectrum, we choose continuum 
wavelengths of 6.8, 7.2, 8.25, 9.1, 10.8, and 11.8 /^m (red points). A spline interpolation was used to determine the level of the continuum 
between these wavelengths (red line). The equivalent width was taken as the integrated line flux above this continuum divided by the 
flux density of the continuum at the line center (blue points). Panels (b) though (d) show selected comparison samples (black line) with 
the median bump DOG spectrum overplotted as a dotted line. The comparison samples have composite spectra that are very similar to 
the bump DOG composite. However, the SMG spectrum in panel (c) shows a stronger power-law component, which could represent dust 
heated either by star formation or AGN activity. 
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Fig. 9. — Illustration of selection effects. Panels (a) and (b): Expected IRAC colors for three local template sources if they were placed 
at various redshifts. Panel (c): The expected 24 fim flux density of three local template sources as a function of redshift. The shaded 
histogram shows the redshift distribution of the b ump DOGs presented i n this paper. The open histogram shows the redshift distribution 
of the power-law DOGs from Houck et al. (2005) and Wccdrnan ct al. (2006). While the IRAC selection allows for a greater range of 
redshifts than we observe, the MIPS bump DOG cut selects strongly for sources at z fa 1.9. This is because bump DOGs tend to have 
strong PAH features, which enter the 24 /im bandpass at 2 S3 1.9. The redshift distribution of the power-law sources is more broad because 
these sources do not have strong PAH features. 
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Fig, 10. — PAH equivalent width versus rest-frame 24 fim luminosity. The bla ck points represent the 107 ULIRGs from [Desai et aEI 
ll2007aD . The green squares represent the average PAH-rich ULIRG spectrum from ISaiina et al.l 1 120071 ) , and the red circles represent the 
median spectrum of the z 2 bump DOGs presented in Figure [7] The cndpoints show the minimum and maximum luminosities of the 
sources that went into the composite spectrum. All equivalent widths were computed using the method illustrated in Figure [8] Compared 
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a The local template spectrum that provided the best fit to the 
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573±52 


<84 


5 


26.1±0.3 


24.7±0.3 


24.4±0.2 


21.2±0.4 


19 


2±0 


1 




43.1±2.7 


52.3±3.3 


45.7±16.6 


33.2±14.8 


529±37 


132±28 


6 


>26.4 


25.0±0.3 


24.3±0.2 


21.4±0.2 


18 


9±0 


1 


>18.0 


40.8±2.7 


56.9±3.4 


67.6±16.7 


15.1±1 1.8 


764±58 


123±28 


7 


25.8±0.2 


24.5±0.3 


23.6±0.1 


20.5±0.2 


19 


0±0 


1 


>18.1 


37.1±2.6 


48.5±3.3 


60.1±16.7 


36.6±14.8 


524±35 


<84 


8 


>26.7 


>24.9 


24.2±0.2 


22.5±0.5 


19 


5±0 


1 


>18.1 


30.4±2.6 


38.0±3.2 


49.6±16.6 


42.0±14.8 


586±40 


300±28 


9 


25.0±0.1 


24.8±0.3 


23.5±0.1 


21.0±0.2 


18 


8±0 


1 


> 17.6 


45.4±2.7 


58.5±3.4 


70.0±16.7 


59.4±14.8 


603±36 


195±28 


10 


25.8±0.2 


24.3±0.2 


23.5±0.1 










> 17.9 


31.9±2.6 


39.9±3.2 


49.8±16.6 


49.2±14.8 


556±41 


<84 


11 


25.5±0.1 


24.5±0.2 


24.3±0.3 










>18.2 


56.8±2.8 


63.7±3.4 


60.2±16.7 


56.9±14.8 


502±37 


265±28 


12 


25.1±0.1 


24.6±0.3 


23.1±0.1 


21.2±0.3 


19 


4±0 


1 


> 18.2 


60.0±2.8 


71.3±3.5 


94.0±16.9 


69.3±14.9 


711±35 


<84 


13 


26.9±0.2 


25.2±0.3 


23.8±0.1 










>17.9 


51.1±2.8 


65.1±3.4 


78.2±16.8 


48.1±14.8 


524±48 


<84 


14 


>25.4 


26.3±2.7 


23.5±0.2 










>18.3 


43.3±2.7 


48.9±3.3 


75.5±16.8 


50.7±14.8 


507±47 


<84 


15 


26.0±0.2 


24.7±0.1 


23.8±0.1 










> 17.8 


30.2±2.5 


38.8±3.2 


53.1±16.6 


44.7±14.8 


601±48 


<84 


16 


>26.6 


26.0±0.5 


24.3±0.2 












57.7±2.8 


94.0±3.7 


182.8±17.6 


262.0±15.5 


535±48 


<84 


17 


>26.6 


26.8±1.1 


24.0±0.3 


21.2±0.3 


19 


2±0 


1 




38.2±2.6 


44.5±3.3 


39.9±16.5 


53.5±14.8 


540±49 


<84 


18 


25.9±0.2 


24.6±0.1 


23.7±0.1 












24.2±2.5 


24.0±3.1 


43.6±16.6 


34.0±14.8 


530±36 


<84 


19 


25.9±0.2 


24.7±0.3 


24.2±0.2 


21.9±0.4 


19 


3±0 


1 




36.8±2.6 


36.0±3.2 


42.7±16.5 


58.5±14.8 


502±44 


<84 


20 


>26.7 


26.4±0.9 


24.3±0.3 


21.1±0.3 


19 


3±0 


1 


>18.3 


16.3±2.4 


22.1±3.1 


33.0±16.5 


30.8±14.8 


501±43 


<84 


21 


>26.5 


27.7±2.1 


>24.7 












21.4±2.5 


25.8±3.1 


10.1±16.5 


27.3±14.7 


508±60 


259±28 


22 


26.4±0.2 


25.3±0.4 


>24.6 










>18.3 


24.2±2.5 


24.7±3.1 


34.4±16.5 


20.4±14.7 


598±50 


<84 


23 


>26.0 


26.7±2.4 


24.7±0.3 












2.9±2.3 


13.0±3.0 


32.5±16.5 


81.3±14.9 


636±49 


<84 



u 



NO TE. — Source 1 through 15 were also the targets of pointed Spitzer/MIPS observations, all of which led to 3tr upper limits of 5000 fiJy and 40,800 /xJy at 70 and 160 /xm, respectively. 
Sec l2~4l 

The ellipses indicate that the source was not covered by the observations. 
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TABLE 3 
Luminosities 



Num 


,,T (0Atim\ 

10 L Q 


Ls— lOOOum (IRS) 

1 nl2T 
10 L 


^8— lOOOum ^raaio,:5ivi^ 

i nl2T 
10 Lq 


Ls— lOOOum (rad: 
10 Lq 


1 


0.98 


7.64 


<3.01 


<11.44 


2 


1.30 


10.17 


<3.42 


<13.01 


3 


0.88 


6.86 


<2.66 


<10.12 


1 


0.87 


6.79 


<2.23 


<8.48 


5 


0.83 


6.44 


3.91 


14.87 


6 


1.39 


10.82 


4.21 


15.99 


7 


0.81 


6.29 


<2.44 


<9.27 


8 


0.97 


7.58 


9.50 


36.13 


9 


0.92 


7.21 


5.09 


19.34 


10 


1.22 


9.48 


<3.32 


<12.63 


11 


1.10 


8.56 


10.48 


39.84 


12 










13 


1.04 


8.13 


<1.53 


<5.81 


14 


0.81 


6.31 


<2.07 


<7.87 


15 


1.01 


7.83 


<2.68 


<10.18 


16 










17 


0.86 


6.68 


<2.54 


<9.65 


18 


0.82 


6.39 


<2.45 


<9.33 


19 


0.77 


5.96 


<2.21 


<8.42 


20 


1.27 


9.87 


<3.66 


<13.92 


21 


6.79 


52.92 


23.85 


90.66 


22 


1.00 


7.80 


<2.68 


<10.18 


23 











a These values were calculated from the observed radio fluxes assuming the FIR-radio correlation calibrated on 
SMGs b These values were calculated from the observed radio fluxes assuming the FIR-radio correlation calibrated on 
nearby galaxies. 



